THz radiation by the frequency down-shift of Nd:YAG lasers 
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The interaction between an intense laser and a relativistic dense electron beam propagating in the 
same direction could down-shift the laser frequency. This process, which can be used to generate a 
coherent THz radiation, is theoretically analyzed. With a set of practically relevant parameters, it 
is suggested that the radiation energy could reach the order of 1 mj per shot in the duration of 100 
pico-second, or the temporal radiation power of 10 MW. 
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A THz light has beneficial characteristics for a range 
of applications (H-Q, as it passes through various non- 
conducting materials while it remains non-ionizing and 
non-invasive because of its relatively low photon energy. 
Most of the recent research efforts have been focused on 
the light source of 1 to 10 THz frequency range 
However, the currently available technologies 
have limitations in generating an intense and short THz 
light source [5[- On the other hand, great advances have 
been made in developing ultra-short and high power en- 
ergy sources, in the context of the inertial confinement fu- 
sion [l2|-[l7j ■ The density of the relativistic electron beam 
can be achieved up to 10 cm , and a 1 /im-light 

source of the intensity exceeding 10 18 W/cm 2 with the 
duration of pico-seconds are currently available. We note 
that an intense, coherent THz light source could be gen- 
erated if the frequency of the above-mentioned lasers gets 
properly shifted to the THz regime. There have been suc- 
cessful attempts to generate the THz radiation using the 
laser driven electron beams [20T - |2^ . 

In this paper, one frequency down-shifting scheme is 
theoretically investigated, which is inspired by the prin- 
ciple of the free electron laser (FEL) utilizing an intense 
electron beam and an appropriate laser. In the con- 
ventional laser-based FEL, a well-collimated relativistic 
electron beam propagates in the opposite direction to 
the intense laser. The amplified wave of the frequency 
lu = 27 2 cjo gets radiated in the same direction with the 
electron beam, where ujq is the original laser frequency 
and 70 is the relativistic factor of the beam, a function of 
the beam speed and the speed of light c. On the contrary, 
if both the relativistic electron beam and the laser prop- 
agate in the same direction, the FEL mechanism would 
result in the wave with the down-shifted frequency of 
uj = o;o/47q, propagating in the opposite direction. Our 
analysis based on the practically relevant parameters re- 
veals that the THz wave with the energy more than 1 
mJ could be radiated per shot. In the rest of the paper, 
we study the interaction between the co-traveling rela- 



tivistic electron beam and an intense laser, estimate the 
amplification efficiency, and then discuss the attractive 
features of this scheme. 

Let us consider an intense laser and a relativistic elec- 
tron beam, both traveling in the positive z-direction. 
We assume the laser is linearly polarized, given by E = 
Eq cos(fcoz — wq£)x and B = Eq sin(fcoz — ujQt)y. In this 
paper, we assume that wo/fco = c so that |E| = |B|, 
which is valid when the electron beam density is much 
lower than the critical density of the laser. This is the 
case for our regime of interest. For simplicity, we ignore 
the electron motion in y-direction. The motion of a rel- 
ativistic electron is given as 
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where y(v)~ 1 = yl — (v/c) 2 and v = |v|. The first order 
perturbed motion of an electron with the initial velocity 
v = vqz is given as 
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- (vq/c) 2 and the approximations of 
-i>o/c)|E| and kv^t — uot = k(vo —c)t 



where 7 1 = 
|E + v/cx B| = fl- 
are used. In driving Eq. @, it is assumed that v£~ /c < 1. 
For an intense laser, there is a time dependent velocity 
(and therefore 7) along the direction of the propagation 
(the z-direction); the electron experiences the eight fig- 
ure motion in the x and z direction. For t4 /c < 1, this 
motion will be time-averaged out so that it can be ig- 
nored. However, the general treatment for the case when 

Vx^ /c = 1 is beyond the scope of this paper. 

Now, consider a counter-propagating electromagnetic 
wave given by E = SEt sin(kTZ + uixfyx and B = 
SEt sin(fcxz + WTt)y. We also assume wr/^r = c. Since 
uj t < wo, the dispersion relationship of the THz wave 
could be significantly different from the one in the vac- 
uum for a high electron density beam. This diviation 
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of the vacuum dispersion relationship would change the 
prediction of the resonant frequency. However, the qual- 
itative conclusion of our analysis does not depend on 
the detailed dispersion relationship of the THz radiation, 
which will be discussed further in Eq. ([7]). The second 
order electron motion is described by 
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to J \mu} cj 7o( w o/ c ) 2 + 7o 

sin(fcoz — u>ot) sin(fc<rz + wji), (3) 

where it is assumed that <C vq so that d,(^v x )/dt = 
(j 3 /c 2 )vl(dv z /dt)+j(dv z /dt). Considering only the res- 
onance term, Eq. ([3]) can be simplified to the form of the 
FEL pendulum equation 

.(2) 



dt \ m c 

cos (4> Q + (k + k T )z - (uj Q - uj T )t) 



(4) 



where </>o is the initial phase, and k(Eq, ko,Vo) is given as 

eE \ 1 



K(E ,k ,v ) 
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The resonance condition of the pendulum equation, 
Eq. (|4]), is given to be 



k T = k (c - v )/(c + v ) = fc / 4 7o- 



(6) 



If the THz dispersion relationship is given as wt(^t) 2 = 
47rn(,e 2 /m e 7o + (?k\ deviating from the vacuum, where 
rib is the electron beam density, the resonant condition is 



inribe 2 /m e is the Langmuir frequency of the beam den- 
sity rib, a- n d f e is the electron distribution function with 
the normalization condition J f e dv z — 1. The Landau 
damping theory suits better than the conventional lin- 
ear FEL analysis [To[ when the electron beam has high 
energy spread, as it is the case here [H, EH, [H, [2(| . 

Eq. © describes the FEL amplification of the down- 
shifted electromagnetic wave. The physical mechanism 
can be understood better by considering the electron 
motion in the reference frame where the electron beam 
is nearly stationary. In this reference frame, the laser 
is seen propagating in the positive z-direction, with the 
down-shifted frequency of ujb = wo/270- If the laser du- 
ration in the laboratory frame is r, the duration in this 
frame is 2jqt. The electron would experience the quiver 
motion by the laser, which results in the electromagnetic 
wave propagating in the opposite direction (negative z- 
direction) unstable. This unstable mode would be am- 
plified and emitted in the negative z-direction with the 
frequency of u> = Wf,/27o and the duration of 2jqt. Seen 
in the laboratory frame, it propagates in the negative 
z-direction with the duration of 47qT and the further 
shifted-down frequency of oj = ujb/Zjo — ^o/47o- For 
instance, if the wavelength of the original laser is 1 /im 
and the duration is 1 pico-second, the wavelength and 
the duration of the amplified wave would be 100 /im and 
100 pico-second for 70 = 5. 

Let us estimate the efficiency of our FEL scheme. 
As an example, consider l-fim Nd:YAG laser with 
the intensity I and the duration r. Denoting lis — 
7/(10 18 W/cm 2 ) and the energy spread of the electron 
beam as SE/E = (, a crude approximation leads to 
(df/dv) u/q = (l/c 2 )(7o/C 2 )- As k in Eq. © can be es- 
timated to be k 2 = (l/107o)ii8, Eq. @ can be re-casted 
into 
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where w^ pc = (4irnbe 2 /m) and we need solve for for 
kx- The exact resonant condition will be different from 
Eq. ([6]), but the down-shift will be the same order as in 
the case when ckr = ujt- As mentioned, the qualitative 
conclusion will not change in any significant way. 

Equating the energy dissipation rate of the electrons 
to the wave growth rate, 
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leads to the Landau damping instability rate (i.e., the 
FEL growth rate) 
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where () is the ensemble average over the electron distri- 
bution, q = ko + hr — ko, and u) = ujq — cut — ojq, <^ pe = 
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where r_i2 = r/(10~ 12 sec), and n%8 — rib/(10 18 cm -3 ). 
For a simple estimation of the plausible physical parame- 
ter, let us assume that the electron beam duration is the 
same as that of the laser and that r_ 12 = 1, C = 0-01 
and 70 = 5.0. Then, the condition 7V > 1 is given as 
^I8 n i8 > 0.003, and the radiation frequency is 3 THz. 

Consider two examples to highlight two competing as- 
pects of this scheme. In the first case, assume lis — 
0.001, nis = 10, 70 = 5.0 and C = 0.01, which satisfies 
7V > 1. Assuming the spot size of both the electron 
beam and the laser are roughly 0.01 cm, the total energy 
of the electron and laser beam is 0.1 J and 0.1 J, respec- 
tively. The possible total radiation of the THz radiation 
per shot would be roughly 0.001 J, assuming 1% of the 
beam energy (£ = 0.01) is converted into the THz radia- 
tion. The radiation time is 100 pico-second and the fre- 
quency is 3 THz. In the second example, assume I\s — 1, 
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ni8 = 0.01, 70 = 5.0 and £ = 0.01, which also satisfies 
the condition 7V > 1. For the same spot size and the 
duration as in the first case, the total energy of the elec- 
tron beam is 0.1 J and the total energy of the laser beam 
is 100 J; 1 ytiJ can be radiated in the THz frequency, as- 
suming the same conversion rate of 1%. The radiation 
time is 100 pico-second and the frequency is 3 THz. 

Two examples given above illustrates the following. 
For a denser electron beam, the total radiated energy is 
higher and the requirement for the laser intensity is less 
severe. In general, a denser electron beam is preferred. 
However, a denser beam would expand very fast via the 
space-charge effect. In order to avoid this, the beam 
needs to be designed to propagate through a plasma of 
comparable electron density to the beam. However, if the 
background plasma density is too high, exceeding the or- 
der of 10 17 cm~ 3 , the resulting THz radiation would not 
propagate any more. It gets reflected toward the inside 
the beam and would not propagate through the beam 
because the THz radiation cannot (can) penetrate the 
beam propagating in the same (opposite) direction. This 
trapping mechanism reduces the total energy radiated 
in the THz frequency. The optimal parameter regime 
would be such that the trapping of the THz radiation 
by the background plasma is not severe, while the elec- 
tron beam is still dense enough for mitigated intensity 
requirement of the laser. Our rough estimate suggests 
that the energy of the laser and the electron beam is a 



few J and the electron beam density (comparable to that 
of the background plasma) is n e = 10 18 ~ 19 cm~ 3 . A more 
detailed calculation of the optimal regime is the topic of 
our future research. 

To summarize, a scheme for the coherent THz radia- 
tion, based on the frequency down-shift of the FEL where 
both the electron beam and the laser propagate in the 
same direction, is theoretically investigated and its prac- 
ticality is examined. In this scheme, the laser frequency 
would be down-shifted by a factor of 4<y 2 , and the total 
radiated power is between 1 fiJ and 1 mJ, depending on 
the laser intensity and the electron beam density. The 
temporal power may reach as high as 100 MW. 

The scheme proposed has many disadvantages. Firstly, 
the required laser intensity is high. Secondly, the required 
electron beam density is also high. The space charge ef- 
fect could be detrimental to the maintenance of the beam 
quality for sufficiently long time. Thirdly, the total ra- 
diated power is not as large as in the case of the gy- 
rotron. Lastly, there are many physical effects that we do 
not consider, which could prevent the scheme from being 
practical. The time dependent electron motion in the z- 
direction and the modified dispersion relationship of the 
radiated THz wave are among those physics. However, 
as the laser technologies develop in a very fast phase, the 
current scheme could get more practical since the major 
obstacles of the conventional FEL, such as the expensive 
magnets, can be avoided. 
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